Lattice approach to the AI = 1/2 ruleQ 
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ABSTRACT 



We discuss a number of old and new methods for computing K — > tttt amplitudes 
on the lattice. They all involve a non-perturbative determination of matching 
^ . coefficients. We show how problems related to operator mixing can be greatly 

l/"") \ reduced by using point-split hadronic currents. 
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g ■ 1 Introduction 

o , 

One of the major puzzles remaining in hadronic physics is the so-called "A/ = 1/2 rule" 
in non-leptonic kaon decays. Decays in which isospin changes by AI = 1/2 are greatly 
enhanced over those with AI = 3/2. For instance, one finds experimentally 
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Although the origin of this large enhancement is not theoretically well understood, we do 
know that, in a QCD-based explanation, most of the enhancement must come from long 
distance, non-perturbative physics. 

Let us briefly discuss the source of the difficulties in calculating A(K — > tttt) in lattice 
QCD i - !]. 

For scales below M w , but above the charm quark mass, the AS = 1 part of the effective 
weak Hamiltonian can be written as 



G " C + (ji,M w )0M(jM)+C4fM,M w )0^(ji) , (2) 



o (±) = [(*tf<0(*yf«) ± - [« <- c l . ( 3 ) 

where 7^ = 7^(1 — 7s)/2, \ u = V u dV* s , Gf is the Fermi constant and \i is the subtraction 
point. 

The operators O^' have different transformation properties under isospin. In particular, 
O^' is pure / = 1/2, whereas OW contains parts having both 7=1/2 and I = 3/2. An 
explanation of the AI = 1/2 rule thus requires that the K — > tttc matrix element of C-O^ 
be substantially enhanced compared to that of C + 0^ + \ 

1 Talk given at "Cortona 1998 - Convegno informale di Fisica Teorica" (Cortona, Italy). 



Part of the enhancement is provided by the ratio of Wilson coefficients, C_/ C+, in their renor- 
malization group evolution from /i ~ M w down to /i ~ 2 GeV, where one finds \C-/C+\ ~ 2. 
This factor is, however, too small by an order of magnitude to explain the AJ = 1/2 rule. 
The remainder of the enhancement must come from the matrix elements of the operators, and 
these are the quantities that we wish to evaluate on the lattice. Attempts in this direction 
date back to the works of ref. |L| 

There are two major problems which arise in the calculation of O^' matrix elements in 
lattice QCD. 

1. Decay amplitudes into two or more particles cannot be directly acccessed in Euclidean 
space, as a consequence of the Maiani and Testa (MT) no-go theorem ||, except in 
kinematical configurations where final state interactions are absent. 

2. Operators mix with operators of lower dimension with coefficients which diverge 
as inverse powers of the lattice spacing. These contributions must be computed non- 
perturbatively and subtracted, leading to prohibitively large statistical errors 0. 

In this note we briefly revisit the lattice approach for the case of Wilson fermions (either 
improved or not) |, recalling old methods and reviewing new proposals ||. 



2 Old and new approaches 

There is a long list of methods proposed in the literature to deal with the problem of 
computing the hadronic matrix elements that enter in the non-leptonic weak decay amplitude 
on the lattice. They are all aimed at bypassing the two major difficulties illustrated in 
the previous section. Here we briefly highlight the merits and the drawbacks of the most 
promising among them. 

1) In ref. H it was suggested to work with m s = and calculate the Euclidean amplitude 



4tU, = (7r(p 1 = 0)7r(p 2 = 0)|OW(/i)|^(p K = 0)) 



(4) 



with all three particles at rest. Using GIM and CPS-symmetry (CPS= CPx s <-> d inter- 
change 0), it can be shown that setting m s = causes all mixings to vanish, removing 
the need for subtractions. Working with the two pions at rest solves the problem of final 
state interactions, as they vanish at threshold. The method requires, however, a large ex- 
trapolation from the unphysical point m s = rrid to the physical one with the use of chiral 
perturbation theory [T^ . 



2) An alternative method || consists in working with a non-perturbatively 0(a) improved 



fermion action (for which there are no errors of 0(a) in the spectrum |] ij and on-shell matrix 
elements of improved currents obey the continuum chiral Ward identities up to 0(a 2 ) [[12] )• 
Choosing quark masses such that rnx = Zm-K, one measures the K — > ttit amplitudes with all 
particles at rest. Since with this choice of masses the momentum transfer vanishes (Ap = 0), 
one can prove that the matrix element of the dangerous, divergent, subtraction is now of 
0(a) rather than of 0(1) and vanishes in the continuum limit. Furthermore pions are at rest 

2 For the use of staggered fermions see ref. 0. 



and the no-go MT theorem does not apply. Besides the need of a (perhaps less severe) chiral 
extrapolation, a problem here is the difficulty of tuning with a sufficiently good accuracy the 
quark masses to have the condition tuk = Ivn^ satisfied. 

3) The old proposal of ref. makes use of the K — > tt matrix elements of (the positive parity 
part of) the weak Hamiltonian and relies on chiral perturbation theory, in the form of Soft 
Pion Theorems (SPT's), to connect them to K — > tttt amplitudes. Since only single-particle 
states are involved, there are no problems with the MT theorem. The disadvantage of the 
method is that the operator mixing problem is much more severe than for the negative parity 
part of TC^g (which enters in the K — > tttt amplitudes), thus making an accurate evaluation 
of the matrix elements of the renormalized operator very difficult. 

The relation between K — > tttt and K — > tt amplitudes is easily found exploiting SPT's. At 
leading order in chiral perturbation theory the physical amplitude takes, in fact, the form 
(for Ap = 0) 
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As the coefficients 7^ ± - ) appear also in the expression for the K — > tt matrix element 
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(ir + (p)\OM(ri\K+(q)) = -tfW -f + 7 (± V q. (6) 

J IT 

by studying this matrix element on the lattice as a function of p ■ q, one can, in principle, 
determine 7^, from which we obtain the K — > tttt amplitudes. 

In order to construct the finite renormalized lattice operators it was suggested in ref- 
erence to use perturbation theory to determine finite mixing coefficients and subtract 
non-perturbatively operators with divergent mixing coefficients. In the case of the parity 
conserving part of O^' that contribute to eq. (§), this can be elegantly done by adjusting 
the only divergent coefficient until the momentum independent part of the K — > tt matrix 
element (8^' in eq. (§)) vanishes. 

4) A method which, in principle, avoids all the difficulties caused by mixing with lower 
dimension operators and which automatically gives the effective weak Hamiltonian with the 
correct normalization is based on the idea of studying at short distances (a |se| ^-qcd) 
the x-behaviour of the OPE expansion of two hadronic currents on the lattice. 

The method does not use chiral perturbation theory, and thus applies equally well to the 
AS = 1, AC = 1 and AB = 1 parts of the weak Hamiltonian. In addition, it allows one to 
construct an improved weak Hamiltonian (i.e. one having errors of 0(a 2 )), if the improved 
version of the weak hadronic currents (1121 are used. 



The method is speculative in the sense that it is likely to require more computational power 
than is presently available, although we expect it to become practical with the advent of 
Teraflop machines. 

We recall that the standard construction of the non-leptonic weak Hamiltonian begins with 
the formula 

«5 = 9 2 w j d±xD w (x;M w )T(j pL (x)Jl L (0)) , (7) 

where 



D w (x;M w ) = J d A p 



P 2 + M 2 



is the longitudinal part of the W-bosoia propagator and J p l is the (left-handed) hadronic 
weak current. One then introduces the operator product expansion (OPE) in the r.h.s of 
ec L- (0)? which is justified by the observation that the dominant contribution to the integral 
comes from distances |x| M' 1 . For physical amplitudes, one obtains in this way 

(h\H%\ti) = ^Y,^,M w )M^(h\0^^)\h'} , (9) 

where di is the dimension of the operator O^(fi), and the functions Cj(/i, M w ) result from 
the integration of the Wilson expansion coefficients, Cj(x; /i) (defined in eq. ( |TT]) below), with 
the TV-propagator, 

a(/i,MjM^ = J d 4 xD w (x-M w )c t (x-fi) . (10) 

The 0^(/j) are quark and/or gluon operators renormalized at the subtraction point /i. The 
functions Cj(/i, M w ) are evaluated in perturbation theory and their running with /i is dictated 
by the renormalization group equation which follows from the /^-independence of the l.h.s. 
ofeq. (0). 

The sum in the expansion (|9]) is over operators of increasing dimension. As the operator 
dimension of Ti^g is 6, we will have to consider in the following only operators with dimensions 
di < 6, since the contribution from operators with di > 6 is suppressed by powers of 1/M W . 

All the intricacies of operator mixing in the definition of the finite and renormalized op- 
erators, OW((u), come about because the integrals in ([7]) and ( [10]) are extended down to 
the region of extremely small x. The complicated mixing of the 0™(^)'s in terms of bare 
operators arises from contact terms when the separation of the two currents goes to zero, 
i.e. when \x\ is of the order of a (the problem is particularly bad here, because chiral sym- 
metry is broken by the lattice regularization). This observation suggests that a simple way 
to circumvent these difficulties is to directly determine the matrix elements of renormalized 
operators by enforcing the OPE on the lattice for distances \x\ much larger than the lattice 
spacing a, but much smaller than Aq^ d , i.e. in a region where perturbation theory can be 
used to determine the expected form of the OPE. 

We imagine proceeding in the following way. If J p l is the appropriately renormalized (and 
possibly improved) finite lattice current operator, we measure in Monte Carlo simulations 
the hadronic matrix element < h\T(J pL (x)Jp L (0))\h'), as a function of x in the region a <C 
\x\ <C ^qcd- The numbers (hlO^dxjlh') entering in eq. (Q) are extracted by fitting the 
x-behaviour of {h\T(J P L(x)Jp L (0))\h') to the OPE formula 

(/i|T(j pi (x)J^(0))|/ i / ) = E^(^^)( /i l C>(i) (M)l^ , ) ) (11) 

i 

where the Wilson coefficients Ci(x; /j,) are determined by continuum perturbation theory using 
any renormalization scheme we like. The scale [i should be chosen so that l//i too lies in 
the range a <C l//i ^qcd- Since we only consider operators of dimension 6 or lower, the 
lattice T-product differs from the right-hand side of eq. ( |TTD by terms of 0(|o;| 2 Aq CD ), which 
is then an estimate of the size of the systematic errors intrinsic in this procedure. As a last 
step we insert the numbers {h\0^ (fx)\h') , determined in this way, in (^), finally obtaining 
an explicit expression for the matrix elements of 7~C^- 



The procedure illustrated above requires the existence of a window, a <C \x\ <C A.qq D , in 
which the distance between the two currents is sufficiently small that perturbation theory 
can be used, but large enough that lattice artifacts, which are suppressed by powers of a/x, 
are tiny. For such a window to exist we need to have an adequately small lattice spacing. 
At the same time the physical volume of the lattice must be sufficiently large to allow the 
formation of hadrons. 

A few remarks may be useful at this point. 

• The method determines directly the "physical" matrix elements of the operators ap- 
pearing in the OPE of the two currents, i.e. the matrix elements of the finite, renor- 
malized operators 0^(/i), without any reference to the magnitude of the PU-mass. 
This means that it will not be necessary to probe distances of 0(1/M W ) with lattice 
calculations. 

• Since it is the continuum OPE which determines the operators appearing in the lattice 
OPE (|^), these are restricted by the continuum symmetries. This is because, for 
\x\ ^> a, the lattice OPE matches that of the continuum with discretization errors 
suppressed by powers of a/x. 

• The /i-dependence of the matrix elements of the operators O^(fi) is trivially induced 
by that of the (perturbative) Wilson coefficients, Cj(x;/i). It compensates the related 
/i-dependence of the functions Cj(/i, M w ) in such a way that the l.h.s of eq. (0) is 
independent of the choice of the subtraction point. 

• Unlike the methods discussed before, this approach automatically yields hadronic am- 
plitudes that are properly normalized (in the renormalization scheme in which the 
Wilson coefficients appearing in eq. ( |TT1) are computed). 

As for the applicability of this method, we remark that, fortunately, in the case at hand there 
are no operators of dimension lower than 6. If lower dimension operators were present they 
would dominate at short distances, since their Wilson coefficients would diverge as powers 
of 1/x (up to logarithmic corrections). In this situation it would be virtually impossible to 
pick out the matrix elements of the interesting dimension 6 operators. 

Operators of dimension 6 have, instead, Wilson coefficients which vary logarithmically with 
x. At leading order their expression is 

(*) 
Zo_ 

/ a s (l x) \ 2 "o / 10 x 
Ci(x;n) oc I ~^-hr 1 =1 + ^o WJ + ..., (12) 

where 7q is the one-loop anomalous dimension of the operator and [3q is coefficient of 
the one-loop term in the /3-function. 

In the case of the AS = 1 part of Ti^, the operators which can appear in (|9|) are 
(eq. (0)) and in addition 

O' = {ml-ml)7s{X-%Kd. (13) 

The GIM mechanism causes O' to vanish when m c = m u , while chiral symmetry requires 
both the quarks to be left-handed and that the GIM factor be quadratic in the quark masses. 
Since O' has dimension 6, its coefficient function depends only logarithmically on x. 



The anomalous dimensions of the three relevant operators (O^ and O') are 

7 J +) =4, 76^ = "8, To = 16- (14) 

Actually, the contribution of O' to the r.h.s. of eq. (|IID can be determined separately, 
since its matrix element does not require any subtraction, and needs not be fitted. The 
anomalous dimensions of the operators O^' are well separated from one another, so it 
might be possible to determine the amplitudes (h\0^ (fi)\h') and obtain the physical matrix 
elements of H^ ff s=1 . 



3 Conclusions 

In this note we have reviewed a number of new approaches aimed at studying the Al = 
1/2 rule on the lattice, using Wilson-like fermions (similar methods could also be used for 
staggered fermions). 

We have reevaluated the methods of refs. |J and || involving K — > titt and K — > n ampli- 
tudes, respectively. The last approach is likely to be more difficult than the first one, because 
of the large number of mixing coefficients which have to be determined non-perturbatively. 
It may however give complementary information to the results obtained with the K —>■ tttt 
methods and provide a check of the accuracy of chiral relations. 

The recently proposed approach || based on the study of the short distance behaviour 
of the OPE of two hadronic weak currents on the lattice is theoretically very appealing. 
The advent of Teraflop Supercomputers may render it a viable method to directly extract 
physical amplitudes from Monte Carlo data. An interesting feasibility study in this direction 
has been undertaken |T3[ . It consists in the investigation of the small x behaviour of the OPE 



of two currents in the two-dimensional 0(3) a- model. A preliminary analysis of Monte Carlo 
data indicates that the measured small x behaviour of the one-particle matrix elements of 
the product of two currents matches the logarithmic behaviour expected from perturbative 
calculations, thus allowing the (non-perturbative) evaluation of the matrix elements of the 
operators appearing in the OPE. 
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